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The magnetic susceptibility of systems from a class of integrable models for doped 
spin-5* Heisenberg chains is calculated in the limit of vanishing magnetic field. For 
small concentrations Xh of the mobile spin- (S — 1/2) charge carriers we find an ex- 
plicit expression for the contribution of the gapless mode associated to the magnetic 
degrees of freedom of these holes to the susceptibility which exhibits a singularity 
for Xh — > for sufficiently large S. We prove a sum rule for the contributions of 
the two gapless magnetic modes in the system to the susceptibility which holds 
for arbitrary hole concentration. This sum rule complements the one for the low 
temperature specific heat which has been obtained previously. 
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1 Introduction 



Essential insights into the properties of low-dimensional correlated electron systems have been 
gained based on studies of integrable models. Exact results on the low energy spectrum of the 
one-dimensional Hubbard and t—J models have made possible their identification as microscopic 
realizations of so-called Tomonaga-Luttinger liquids allowing a complete classification of their 
critical exponents Recently, the t—J model - introduced to describe the effect of hole 

doping on an S = 1/2 Heisenberg antiferromagnet - has been generalized to larger values of 
S j|,|5| giving the framework for detailed studies of quantum effects in the double exchange 
model. Fine tuning of the coupling constants in these models describing spin-,!? chains doped 
with mobile spin-S" = (S — 1/2) carriers (called 'holes' below) an exactly solvable model has 
been constructed ||, 0] which interpolates smoothly between the spin-5 1 (spin-S") Takhtajan- 
Babujian chain [PHT0||. 

An important open question is the proper identification of effective field theories describing 
the low energy/low temperature behaviour of these systems ||, 0,^,0. In the limit of low hole 
concentration the low temperature specific heat indicates that the two gapless magnetic modes 
found for vanishing magnetic field can be described in terms of a level-2S' SU(2) Wess-Zumino- 
Witten (WZW) model (which is the effective field theory of the undoped Takhtajan-Babujian 
spin-5 1 chain P~p!0f|) and the minimal model M.2S+I1 respectively ||[7|]. Upon variation of the 
doping a continous transition to a product SU(2) 2 s-i <§) SU(2)i satisfying a sum rule in the 
coefficients of the low temperature specific heat is observed . For S = 1 a field theoretical 
description of these degrees of freedom has been proposed based on four Majorana fermions H] 



(see also [TT],[T2|]). However, to determine the coupling constants present in this continuum 
theory additional insight into the physical properties of the massless excitations, in particular 
their response to an external magnetic field, is necessary. 

In this paper we analyze the Bethe Ansatz equations to compute the zero temperature 
magnetic susceptibility of the integrable doped Heisenberg chains in the limit of vanishing 
magnetic field H — > 0. Below we briefly review the thermodynamic Bethe Ansatz for the model 
to provide the foundation for the subsequent analysis. In Sect. § we concentrate on the limit of 
small magnetic fields and derive a particular sum rule for the contributions Xi an d X2S of the 
two gapless magnetic excitations to the susceptibility, similar to the one for the specific heat 
mentioned above. In Sect. |] we formulate a matrix Riemann-Hilbert problem (RHP) whose 
solution determines \i an d X2S individually. It is shown that the symmetries of this RHP 
immediately imply the sum rule mentioned above. In Sect. |5| we present the main result of this 
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paper which is the explicit calculation of %x an d X2S m the limit of small hole concentration. 
The Hamiltonians considered in this paper are of the form 

L 



W (S) = E{^i + ; Ci}- C 1 - 1 ) 



n=l 



Here and Af^ describe the hopping of the holes and the (antiferromagnetic) exchange 
between sites i and j of the lattice, respectively. Similar Hamiltonians arise as effective spin 
models obtained from the double exchange (DE) model |T3|,|T4}j in the limit of strong ferromag- 



netic Hund's rule coupling between the spins of the itinerant electrons and the local moments 
with spin-S' @,|]. In the basis of the relevant spin multiplets the kinetic part of the Hamiltonian 

reads 

= -V lJ Q s (S t -S J ) (1.2) 

where S, is a spin-S 1 (5") operator describing the particle (hole) on site i. The operator Vij 
permutes the states on sites i and j thereby allowing the holes to move and Qs{x) is a polynomial 
of degree 2S — 1 leading to a hopping amplitude depending on the total spin St — §, §, • • • 2S — 
| on sites i and j. In the DE model these hopping amplitudes favour the formation of a 
ferromagnetically ordered state, this feature is not shared by the integrable chains 0,0]. In terms 
of a microscopic model the electronic hopping has to be dressed by bilinears in the localized 
magnetic moments, similar to bond-charge type interactions considered in generalizations of 
the Hubbard model. 

As a consequence of ST/ (2) invariance the exchange operators can also be written as a 
polynomial in Si ■ Sj, the precise form of these polynomials depends on the local configuration, 
i.e. particle-particle, hole-hole and particle-hole exchange are different. As mentioned above, 
the pure limits are the integrable spin-S and S' Takhtajan-Babujian chains determining the 
two former processes. 



2 Bethe Ansatz for the doped spin-S* chain 

Starting from the ferromagnetically polarized eigenstate of the undoped chain we consider 
excitations obtained by adding Nh holes and - in addition - lowering iVj_ spins. The spectrum 
of these states can be studied using the algebraic Bethe Ansatz. They are parameterized in 
terms of Nh real numbers v a and Nh + iVj_ complex numbers Xj solving the Bethe Ansatz 
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equations (BAE) @,[7| 

A, + iS\ u T-r 1 \j — Ai. + i tt A, — v a — 



L N h +N l N h 



A i - iS J i£j ^ ~ Afe ~ 1 a=l A ^ - "« + 5 

j = l,... ,jV & + iV ; (2.1) 

N h +N, < 

1 = r j , a = l,...,iV fc . 

In the thermodynamic limit L — > oo general solutions of these equations are known to consist 
of real hole rapidities v a and complex n- strings of spin- rapidities X]' k = \f ) + |(n + 1 - 2Jfe), 
= 1, . . . , n with real centers A^ . Considering solutions of ( |2.1| ) built from Nh hole rapidities 
and M n A-strings of length n we can rewrite the BAE in terms of the real variables u a and A^. 
Taking the logarithm of (|2.1|) we obtain 



L6 n , 2S (Af) = 2vr jW + £ £ ~ nm (Af - A< m >) - £ £ n (Af - />, 



m=l fc=l a=l 

oo M n 



n=l j=l 



where # n (x) = 2 arctan(2x/n) and 



@nm(%) — @m+n-l {%) + $m+n-3 + • • • + 0| m _ n |+l (x) , 

2 nm (x) = 6» n+m (x) + 26> n+m „ 2 (x) + . . . + 2#| n _ m | +2 (x) + (1 - 5 nm ) 9\ n - m \ (x) . (2.3) 

From Eqs. ( |2.2| ) the allowed values of the quantum numbers I a and are found to be 
(t nm = min(n, m) - 5 nm /2) 



oo 

I 

/J < 



- M n > ^ < 9^ min (^ 2 5) - E *nmM m + - (JV h - 1) . (2.4) 



2 

m=l m=l 



A given choice of these numbers the solution of ( |2.2|) uniquely determines a particular eigenstate 
of the system with energy 

oo M n , . v 

E = E E { e « } ( A ? } ) + ^} - (m + H ^ - LSH ( 2 - 5 ) 

n=l j=l ^ ' 

and magnetization 

^ ot = LS-^nM n + -iV h . (2.6) 



n=l 
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In Eq. (|2.5| ) H is the external magnetic field, p is the chemical potential for the holes and the 
"bare energies" of the A-strings are en\x) = —2-k (A n> 2s * s) (x) with s(x) = 1/(2 cosh ttx) and 

A nm {x) = —E' nm (x) + S nm 5(x) (2.7) 

((/ * 9) ( x ) denotes a convolution). 

In the thermodynamic limit L — > 00 with Nh/L and M n /L held fixed one can introduce 
densities p(x) for the hole rapidities and <r n (x) for the A-strings of length n. The BAE (|2.2| ) 
turn into linear integral equations for these functions 

d n (x) = (An^S * s) (x) - ^ (Am * 0"m) (%) + (»n * P) (%) , 

m 

p(x) + p(x) = ^ (a n * a n ) (x) . (2.8) 

n 

Here, 2iia n (x) = 9' n (x) = An/ {Ax 2 + n 2 ), and p(x), a n (x) are the densities associated with the 
distribution of vacancies of quantum numbers I a and J. in the intervals (|2T4T), respectively. 
Similarly, the energy (|2.5|) and magnetization ( |2.6|) expressed in terms of these densities read 

00 „ „ 
E/L = J2 dxe^(x)a n (x) -p dx p(x) - H 5 t 2 ot , (2.9) 

n=l J J 

andQ 

- 5 t * ot = S - ^2 n I dxa n (x) + \ I dxp(x) = ^ ton / dxa n (x) . (2.10) 

n=l ^ ^ n^oo J 

Analysis of these equations shows that the ground state configuration of the system at 
temperature T = consists of up to three condensates formed by the hole rapidities and A- 
strings of length 1 and 25 only. These "Fermi seas" are characterized by their end points Q 
for p(x) and Aj for <Tj(x), j = 1,25. In the following we drop the distinction of particle and 
vacancy densities by identifying p(x) with p(x) + p(x) and similar for the densities a n (x) of the 
A-strings. This is possible since 

p(x) for |x| > Q 



Using Mj = for j ^ 1, 25 in Eqs. ( j2.4|) and ( |2.6| ) one can express the magnetization in terms 
of (T25 only 

jS^= r dxa 2S {x) . (2.12) 



1 The second expression is obtained by integrating Eq. (eT8J) over all 



3 Thermodynamic properties in a small magnetic field 



For vanishing magnetic field one can show that Aj — ► oo. Hence, in the limit of a small magnetic 
field H which we are interested in below one has Aj 3> Q allowing to eliminate the <Tj from the 
integral equation for the density p of hole rapidities 0. This results in 

rQ 

p(x) - dyR(x - y)p{y) = (a 2S * s) (x) . (3.1) 

J-Q 

The kernel of this integral equation is R(x) = a 2 * (1 + c^) -1 and the concentration of the holes 
is Xh = Nh/L = J_gdxp(x). In this regime the boundaries of integration or "Fermi points" 
±Q are a function of the hole chemical potential p alone: they are determined through the 
condition k,q(±Q) = for the dressed energy kix) of the excitations associated with the hole 
rapidities which in turn is given as solution of a linear integral equation similar to ( | . 1|) 



K {x)- dy R(x - y)n Q (y) = - {2ixa 2S * s(x) + p] . (3.2) 
J-Q 

Excitations with charge rapidities near ±Q are massless. The effective low energy theory for 
the charge mode has been identified as that of a free boson. The velocity of the charge mode 
can be obtained from the dispersion ( |3.2| ) 

1 dn 



(3.3) 



2np(Q) dx 

The integral equations for the relevant densities of A-strings o\ and a 2 s can be rewritten as 



ai ( x ) = I dy s(x - y)p{y) + / dy K n (x - y)a\{y) + / dy K 12 (x - y)a 2S (y) 

'\y\>M J\y\> A 2s 



(3.4) 



°2s(x) = s(x)+ dy K 21 (x -y)a x {y) + / dy K 22 (x - y)a 2S (y) 

'\y\>M Ay\>^2s 



The kernels of the integral operators are easiest given in terms of their Fourier transforms 

1 / sinh(5' — l)u sinh \uj 



2cosh|u; sinh(S'- \)u 1 sinh |w sinh(S - 1)uj + sinh(S' 



2 

1 



2J u 



(3.5) 



Note that for any given positive x the contributions from the intervals y < —Aj in Eqs. (|3.4j) 
can be neglected in the limit H — > 0. After replacing Jj y | >A by these equations are a system 
of Wiener- Hopf (WH) integral equations. Furthermore, the replacement of the driving terms by 
the large-x asymptotics will not change the behaviour of the solutions of these equations near 
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the Fermi points x ~ Aj ^> 1 which contains all the relevant information for the low energy 
properties of the system. Hence we arrive at the following set of equations for the densities of 
the A-strings valid for i>1 

poo poo 

a x {x) = B e~" N + dy K u (x - y)a x {y) + / dy K 12 {x - y)a 2S {y) , 

JM Ja 2S 

(3.6) 

POO POO 

a 2S (x) = e" wN + / dy K 21 {x - y)cn{y) + / dy K 22 {x - y)a 2S {y) . 

J K\ JA 2S 

with B = J® dx exp(irx)p(x). In the same spirit we obtain a system of WH equations for the 
dressed energies of the magnetic excitations 

POO POO 

e 1 (x) = -2nAe-*W+ dy K n (x - y)e 1 (y) + / dy K 12 {x - y)e 2S (y) , 

</Ai Ja 2S 

(3.7) 

poo poo 

e 2S (x) = -H-2Tre- n ^+ dy K 21 (x - y)e 1 (y) + / dy K 22 (x - y)e 2S (y) . 

1 J A x J A 2S 

with 2ttA = — J gdy exp(irx)K (y). The solution of ( ^.7[ ) determines the values of Aj as a 
function of the external magnetic field H by the condition €j(Aj) = 0. 

Both magnetic modes e\{x) and e 2 s{x) allow for massless excitations near these Fermi points, 
in the limit H — > their velocities are 

^ 25 = lim 251 = tt , Vl = hm -^/t = *~5- 3 - 8 

Their critical properties have been studied using the thermodynamic Bethe Ansatz. In the 
limiting cases of the undoped (completely doped) model where the model reduces to the inte- 
grable spin-S" (spin- ( 5 — 1/2)) Takhtajan-Babujian spin chains, respectively, e 2 s (x) reduces to 
the massless spinon mode of these models which can be described by a SU (2) level-2S' (level- 
(2S — 1)) WZW model. For arbitrary doping the contribution of the magnetic modes to the 
low temperature specific heat has been found to be ||[7| 

C ! zC u + C lS %(*Z+*) (3.9) 



3 \V 2S V X 

where the parameters Cj satisfy the "sum rule" 



45-1 . . 

C2S + C1S2 2STT (3 ' 10) 



independent of the doping, i.e. the values of the parameters A, B in Eqs. ( p.6|) and (|3.7|) . In 
the limiting cases mentioned above c 2 s becomes the conformal central charge of the SU(2) 2 s 
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and the SU(2)2s-i WZW models, respectively. From the other parameter, ci, the low energy 
theory for e% is identified as the minimal unitary model .M2S+1 in the limit of vanishing hole 
concentration (xh — > 0) and a free boson with c\ = 1 and compactification radius i? = l/y/2 as 
Xft — > 1 0. The latter is equivalent to the level-1 SU(2) WZW model. 

The magnetic field dependent part of ground state energy of the system for small H ex- 



pressed in terms of the solutions of (|3.7|) is 



e s = -E = -2n[ — dx e^e^x) + - / da; Ae~™e x {x) = -- X H 2 . (3.11) 
L \V2s Ja 2S vi J Al J 2 

This implies, that the zero field magnetic susceptibility x wm have contributions from both 
massless magnetic modes - very similar to the feature found for the specific heat (|3.9| ) 



X = X2S + Xi (3.12) 

A sum rule for the contributions Xj can be obtained by considering formally the case A = B 
which implies V\ = v 2 s = then the solutions of the integral equations ( |3.6| ) and (|3.7|) are 
related through 



Using this relation with the asymptotic behaviour \im x ^ oc e2s(x) = 2S H in ( |2.12j ) one imme- 
diately finds^ 

S 

V2SX2S + V1X1 = - ■ (3.14) 

7T 

From the known susceptibilities for the limiting cases of the spin-S* and spin-(S'— 1/2) Takhtajan- 
Babujian spin |L0|| chains we expect that (^25X25) decreases continuously from S/ir for vanishing 



hole concentration Xh to (S — 1/2) /n for Xh = 1. This implies ViXi = for Xh = 0. However, 
since v± vanishes in this limit as well, these simple considerations do not rule out a finite or 
even singular contribution Xi to the susceptibility ( |3.12|) . 

4 Calculation of the susceptibility 

To calculate Xj we have to analyze the Wiener-Hopf equations ( |3.7|) for the dressed energies 
in the limit H — > 0. In terms of the functions fj(x) = £j(Aj + x) the WH equations can be 
2 Note that these arguments are easily extended to the case of non-zero temperatures. 
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rewritten as (t = A 2 s — Ai) 

/*00 POO 

h{x) = -2nAe^ +x \+ dy K u {x - y)h{y) + / dy K 12 (x - y - t)f 2S (y) 

Jo Jo 



(4.1) 



f2s(x) 



| poo poo 

-H -2ite-^ +x \+ / dy K 21 (x-y + t)f 1 (y)+ dy K 22 (x - y)f 2S (y) . 
1 Jo Jo 



Again, the Aj have to be determined such that fj(0) = 0. By Fourier transformation of (|4.1|) 
we obtain 



G T {u)F + {uu) -F_{uu) = T(u) 



(4.2) 



Here, F± are the two-component vectors 



F±(u>) = ± 



f™dx e ±iux h{x) 
J°°dx e ±iuJx f 2 (x) 



(4.3) 



of analytical functions of to in the upper (lower) half-planes, respectively. G(oS) is a matrix 
given by 

2 (1-e-H) 



G(w) = U~\u) (/ - K (5) (cu)) C/(o;) , det G{u) = (l + e~ M ) 2 — 



(4.4) 



e -(as-i)M) 

(here U(u) = diag (exp(— iuAi), exp(— iuA 2 s)) and / is the 2x2 unit matrix) and T(uj) is the 



vector 



T{uo) = ttH8{uj) 




1 



4tt 2 / A 



(4.5) 



Now the solution of the equation (4.2) can be given in terms of the one of the regular matrix 
Riemann-Hilbert problem (RHP): 



Z(u>) — > I, uj — > oo, 



Z(u) is analytical for a; ^ I 



(4.6) 



Here and below the subscripts ± denote the limit values of functions from the left (right) of 
the contour. 

If the solution of the RHP ( |4.6|) is found, then Eq. (12) can be written in the form 



zl(u) 



F + {u) 



ZT(co) 



Fju) 



Z T _(u) 



T{u), 



(4.7) 



which is evidently solved by 



1 f°° 
F+(u) = —ZI(lu) I du 



Z T _{u) T{u) 



-i 



3 u — UJ + 

Now we can substitute here the explicit form Q4.5| ) of T(oj). The integration of 5- function is 
trivial. For the second term one can shift the integration contour in ( f4.8| ) to the lower half- 
plane. Then the only singularity of the integrand is simple pole at the point u = —in and we 
obtain 



-l / \ liri r i -i 

Zliuj) Z T (-m) Q(-tir) . (4.9) 



The condition fj(x = 0) = implies —i\m\ ul ^ 00 ojF + {ijj) = 0. Hence the boundaries of 
integration Aj are given as a funtion of the magnetic field H and the parameter A by 



H 

4tt' 



tTi 



-III) 



zi(o) 



-i / o 



(4.10) 



Similarly, the ground state energy Q3.ll ) is now 



-2tt 



V2S Jo 



e- nx f 2S (x) + 



Vl Jo 



-27rO(-z7r) J y _1 F + (i7r 



(4.11) 



(V = diag (vx, V2s) is a diagonal matrix containing the Fermi velocities ( |3.8|) of the magnetic 
modes). Using ( ]4.1U|j in ( O ) we can express F + (iix) through the solution of the RHP 



F + (in) = ^Zl(in) 



zm 



-i / 



(4.12) 



and finally have the magnetic susceptibility in terms of Z as 



X 4tt 



Z_(0) 



-i 



Z_{-m) V- 1 Zl{i7i) 



zm 



-i 



(4.13) 



22 



The symmetries the matrix G(u) allow to deduce some general properties of the solution Z 
of the matrix Riemann-Hilbert problem (|4.6| ). As a consequence of the identity G T (— uj) = G(u) 
one has Z + (oj) Z t (—uj) = I. This allows to rewrite the expression Q4.13 ) in terms of Z + (0) and 
ZAiir) alone 



1 



\ = ^ \ G-\0) Z + (0) Z T + {m)V~' ( Z + (0) 



-l 



-i 



Z T Hm 



(4.14) 



22 
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Eqs. (U) and (p|) yield 



Replacing V by a unit matrix this immediately reproduces the susceptibility sum rule (|3.14|) . 



In Figure [l] we present numerical results from the integration of Eqs. ( j4.1|) for 5 = 1 and S = 3, 
respectively. In the limiting cases Xh — > and Xh — > 1 corresponding to undoped spin-S 1 and 
-(S — 1/2) chains, the velocity v\ vanishes. The singular behaviour of the susceptibility near 
Xh = 1 is a consequence of the contribution xi ~ ^/(^Vi) to the susceptibility. 



5 Susceptibility at low doping 

Interestingly, however, the dependence of the susceptibility on the hole concentration is found 
to depend strongly on the value of S for Xh — > 0: for all values of S the velocity V\ vanishes 
with Xh- In fact, since Xh <C 1 corresponds to small Q in the integral equations (|3.1| ) and (|3.2| ) 
for the densities and dressed energies of the charge excitations, the latter can be solved by 
iteration giving (ip(x) and tp^(x) are the digamma and polygamma functions, respectively) 

12 4, (>M) - V, (2£tl) » 1 ' 

and = ixA/xh oc x\. As discussed above, we expect V\X\ — > in the limit Xh —>■ which 
agrees with numerical findings. The contribution xi of the 'minimal' mode ei(a;) alone to the 
susceptibility, however, is singular for sufficiently large value of S. 

In general, however, we have v± ^ V2s and to find the susceptibility we need to solve matrix 
RHP (4.6). As in the studies of the low temperature thermodynamics |6|, ^] we make use of 



the fact that we have A 1 for small hole concentration Xh- for small H one finds one has 
OcAi < A25 in this regime. Thus, the jump matrix G(uS) ( |4.4| ) contains the large parameter 
t — > +00. This allows to find the solution of the RHP ( f4.6|) asymptotically and thus to determine 
the susceptibility as a function of A. Since the treatment for S = 1 and S > 1 differs in we 
present the main results of this asymptotic analysis separately below. Further technical details 
are presented in the Appendix. 
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5.1 S = l 

To solve the RHP Q4.6|) we factorize its solution into a product of two matrices 

1 

a(w) 

Here the function a(u) is analytical for uo R, a(c<j) — > 1 as u 



Z(u) = $M 



(5.2) 

oo, and solves the regular 



scalar RHP on the real axis: 



It is easy to see that 



q. + {uj) det G{uj) 

a+(u) (l + e- H y 2 for u G 

2ll / ILO 



(5.3) 



(5.4) 



r 2 (j - § ) V 27r e/ 
and a_(co>) = «^ 1 (— cj). In particular, we have a+(0) = 2 and a + {iit) = ir/e. 

From ( |4.6|) we find that the matrix $(u;) in ( |5.2|) is analytical for u £ R, it approaches 
identity at u — > oo and its boundary values on the real axis are related by 

$_(o;) = $ + (a;)G$(a;), w Gl . (5.5) 



The matrix 



G*(u;) 



1 



a_ [OJ) 
2 cosh * 



-iujt 



Atot 



l + e-H 



(5.6) 



/ 



\ 2 cosh | 

can be factorized into the product G$(u) = M+(u)M-(u) with (cr± are Pauli matrices) 



Af+(u;) 



2 cosh f 



M_(u;) 



2 cosh | 



(5.7) 



Obviously M±(±lj) are analytical in the strip < < it. In particular M±(u) can be 
analytically continued into the strip < ±Qw < a for some a G (0, ir) where the off-diagonal 
entries of M± vanish in the limit t — > oo. Following the asymptotic analysis of a similar problem 
in Ref. |L5]| we introduce a deformation of the RHP (|5.5|) which can be solved in terms of a 
absolutely convergent series in exp (— irt). In Appendix [A.l| the asymptotic behaviour of the 
matrix $(u;) in the points u = iir, is computed from this deformed RHP (p and g are functions 



of t defined in Eq. flQ7D ) 

$ + (i7r) 



1 +P 2 (g + 1) 



P 

1 -p 2 



+ o (p 2 ) 



(5.8) 
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, l + 2p-2p 2 2p , 
$ + (0)= \ \\ +o (p 2 . (5.9) 

1 1 - 2p - 2p 2 l-2p 2 1 

Using Eqs. ( |4.10p , ( [4.14| ) we obtain A = 2irp 2 + o (p 2 ) which finally gives the zero field limit of 
the magnetic susceptibility for small doping (or, equivalently, A <C 1) 

X = - (— + -) , (5.10) 



7T \ V 2 Vl 

where 

^-^( io 4 +2c ) + °(- 4) - < 5n > 

This result refines the estimate given in Ref . . In terms of the hole concentration the contri- 
bution xi to the magnetic susceptibility is 

Xl = -4 (log + 0.217927...) (5.12) 

which is in excellent agreement with the numerical data in Figure |](a). 

5.2 S>1 

Similar to the case S = 1 the solution of the corresponding RHP is given by the product of two 
matrices 

*) = »H| ftH ° I, (5.13) 

ft(u>) 

where solve scalar regular RHPs with canonical normalization condition and 

/3i_(u;) = A+HGnH, wgR, 

/3 2 _H = /M«)^rr' (5.14) 
G n (a;) 

Below, we will need only the ratio ^(uj) / (3i(uj) = a^ 5 ^(c<j), which is equal to 

{s)L ^_{2S-i)^s-D + k r(-g)r(-g(25-i)) / 



«W) = — V-§ v 2?r/ ?%\ ^ -— . (5.15) 

+ ^ ' ( 2 <?)^ +1 r2 (-^) V ^e) V } 

As before a+ (cj)a_ (— 10) = 1 and in particular 

af (0) = 2517, af f |) = ^ = (25 - I)" 1 / 2 . (5.16) 
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The asymptotic analysis of the RHP for 3>(a>) is again based on the factorization of the 
corresponding jump-matrix G$(u>) and subsequent deformation of the original jump-contour 
(see Appendix |A.2| for details). The matrices $+(0) and $ + (i7r) needed for the calculation of 
the susceptibility are 

* + (0)H 1 + "'™-" ?p2/2 UlP | + o(A (5.17) 

(1 — ufp 2 /2)r] — uip l — u\p 2 /2 



and 



i _i_ u Ip 2 uip _ u 2 p 2 
I {S_1) 2 5+1 ft 2 |+«(P 2 )- (5-18) 

5-1 5-2 1 2(5+1) 



Here we have used the following notations 



1 (5) ( mk\ . ( irk\ 



In the case S = 2 the entry <J> 2 i+(^ 7r ) should be replaced by 

$ 21+ (ztt) =u lP -7iu 2 p 2 (t-id z \oga { +\z) )+o{p 2 ), S = 2. (5.20) 

\ Z—vkJ 

This difference, however, does not affect the final results for A which is 

A = ^jP-V^P S+1 + o(p S+1 ) . (5.21) 
[o + l)e 

For susceptibility of the spin-S" model at small hole concentration we obtain 

X -S(1Z1 + 1), (5.22) 



7T \ V 2 S Vi 

with 

S 2 



Substituting here u\ from (|5.19 ) and p from ([5.21 ) we finally obtain 



1 / T 2 (S) ( 2S-\ 



As a function of the hole concentration the contribution of the minimal mode e\ to the suscep- 

o 2-S 

tibility is \i °^ x h +1 > which will give a singu 
Fig. 0(b), xi diverges as for small x^. 



2 ^ — J 

tibility is Xi ck x h +1 > which will give a singularity for 5 > 2. For 5 = 3, the case depicted in 
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6 Summary and Conclusions 



In this paper we have studied the magnetic properties of a class of integrable models for doped 
spin-5* Heisenberg chains in the limit of very small magnetic field. Of particular interest was 
the singular behaviour of the susceptibility as a function of the hole concentration for small 
Xh- The response of the two magnetic modes present in these systems to an external field is 
essential for the proper identification of an effective field theoretical description of the system 
which in turn will allow for analytical investigations of the properties of the metallic phases 
of quasi-one dimensional transition metal oxides. The susceptibility due to the 'background 
mode' 62s{x), which interpolates smoothly between the spinons of the spin-5 1 and -(S — 1/2) 
Takhtajan-Babujian models when the hole concentration is varied between Xh = and Xh = 1, 
decreases monotonically from S/n to (S— 1/2) /n as a function of Xh- The contribution Xi of the 
second magnetic mode ti(x) can be singular near Xh = 0, 1 as a consequence of the vanishing 
of the corresponding Fermi velocity. For — > 1, where few spin- S particles propagate in a 
background of the spin- ( 5 — 1/2) holes this singularity may be understood in a similar way 
as for an underscreened Kondo impurity (although the net impurity moment vanishes in the 
ground state). For small doping the effect of the doping on the susceptibility is more subtle. 
Our analysis shows that lim Xh ^o Xi is non-zero only for S > 2. The effect of the (non-universal) 
velocities on the singular behaviour of the susceptibility ( |3.12|) and the specific heat ( |3.9| ) can be 
eliminated by introducing a "Wilson Ratio" for the relative contributions of the two magnetic 
modes 

Rw = ■ (6-1) 

With this definition the contribution from 62s is interpreted as that of the host while associating 
the mode e\ with the impurities which is certainly justified in the limits Xh — > 0, 1. In these 
limits we find 

35(25+1) {■ n 

R w ~ { P+Wis-i) « x »-+ j (6 . 2) 

25+1 Xh ^ 1 



where according to (5.11) and ( |5.24| ) £ vanishes as x^logXh for 5 = 1 and as x^ 1 for S > 1. 



Extending the interpretation of the two magnetic modes to arbitrary hole concentration one 
finds a continuous change of Rw with Xh (see Figure 0): as was to be expected from the low 
temperature specific heat ( p.9|) the ratio defined in ( |6.1|) is not universal. 

Note that our results are the leading contributions to the susceptibility for small magnetic 
field. Due to approximations such as that of Eqs. ( |3.4| ) by ( |3.6| ) additional terms oc 1/ log H 
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to the susceptibility have been neglected. For the undoped systems the precise form of these 
contributions (which should coincide with the ones present in \im Xh ^o X2s) is given in Ref. [TO 



Finally, we mention that the solution of the matrix Riemann-Hilbert problem ( f4.6|) pre- 
sented in this paper opens the possibility to study the critical exponents in the asymptotics of 
correlation functions in the H = phase of the doped Heisenberg chains. In the case of a simple 
Tomonaga-Luttinger liquid the spin part of correlation functions is completely determined by 
the SU(2) symmetry and only the compactification radius of the boson related to the charge 
mode gives rise to anomalous exponents |L|. In the doped spin- 5* chains considered here one 
should expect an additional doping dependence of the critical exponents due to the presence of 
two magnetic modes: in spite of the SU(2) invariance of the system a subtle balance between 
these modes is preserved which is manifest in the sum rules for the low temperature specific 
heat ( |3.10| ) and the zero field magnetic susceptibility ( p. 14| ) . Further investigations are also 
necessary to eludicate the relation of these quantities in the finite field phase and the one with 
vanishing magnetic field. 
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A Asymptotic analysis of the RHP for large t 



A.l S=l 



Using the factorization of the jump matrix G$(u) in Eq. (5J5) we follow Ref. [TjJ by defining 
a new matrix U (ui) as 

$(u/), > a, 

l 7 U')={ <$>{u)M + {uj), 0<%u<a, (A.l) 
^(u/JMrV), -a < 3w < 0, 

where a G (0,7r) (see Fig. D). Then C/(a;) solves the following regular RHP: 

[/ (u) — > /, c<j — > oo, 

U {oS) is analytical for uo ^ T, (A-2) 
C/_(a;) = U+(u)Gu(u), wgT, 

where we have introduced 

^t/M = i . (A.3) 

The contour r = r + ur_ is shown in Fig. |2| (the arrows show positive direction). It is easy to see 
that U (u)) has no cut on the real axis. Thus LL (u) can be analytically continued inside the strip 
| So; | < a. Similarly, U + (uo) is analytical in the remaining domain > a, where it coincides 
with the matrix (more precisely t/+(c<->) = $+(o;) for Su; > a, but U + {uo) = $_(c<j) for 

Quo < —a). 

From the explicit expression (|5.7| ) we have M± — I ~ e~ ai on the contour r^. This implies 



U(u) = I + o (e~ at ) , t -> oo, (A.4) 



for the solution of ( A.2 ) where a is arbitrary from the interval (0,7r). Thus, for i — > oo the 



asymptotic behaviour of the solution to (|5.5[) is 

$H = / + 0(e^), |^|>vr, 

(A. 5) 

$(cu) = / + O (e^ *) , < | So; | = uo < 7r. 

This accuracy, however, is not sufficient to calculate the leading doping dependence of the 
susceptibility. In particular this estimate can not be used to determine $(0) entering Eqs. (|4. 10Q 
and ([4.141) . To obtain the subleading contributions we can use singular integral equations 
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equivalent to the RHP (|A.2p . Note that only $ + (i7r) = U+{vk) and $+(0) are needed for the 
susceptibility. 

From ( |5.6| ) we find detG$(u;) = 1 and G$(u) = G$(—lj), hence we conclude that 

det $(u;) = 1 , $+(w)$ t (-cj) = / . (A.6) 

Hence we have 



$n + (i7r) = U u+ (iir) = C/ 22+ (-i7r), $ 2 i + (ztt) = C/ 21+ (i7r) = -U 12+ (-m) 

$i2+(«7r) = f/i 2+ (i7r) = -U 2 i+{-iir), $ 22 +(wr) = kW(wr) = U 11+ (-iir) 

at = 277 while $(u; = 0) is parameterized by a single parameter ip 



(A.7) 



, cos 09 + sin 09 sin oj . 
$ + (0) =| * * ^ | (A.8) 
cos 09 — sin 09 cos o? 



as a consequence of Eq. (|A.6|) together with jump condition (|5 
Consider now the singular integral equation for U + (uj): 



y+M = J --L/g±W(^w^ (bi (A . 9) 

27u ./r 2 — u; + 



where u; + means that u; is shifted from the integration contour to the left. In components, the 
equations for the entries {7 2 i+(u;) and t/ 22 +(u;) read: 

1 /" [/ 22+ (,2)a + (>)e i2 ' 

U 2 i+{u}) - 



2m J r 2 cosh § z — c<j + ' 



(A.10) 



U 2 2+{u) = 1 



1 f U 2 i+(z)e- izt dz 



2iri Jy_ 2a_(z) cosh ^ z — co + 
Let w G T_ in the first of ( [A. 10 ). Then shifting the integration contour into the upper half- 



plane, we obtain 

U22+ (in (2k + l))a + (m(2k + l)W t7r ( 2fc+1 ) 

Un+(u) = % V (-l) fc ^ '- h, ; 

21+v J ) u-in(2k + l) 



k=0 



(All) 



Since t/ 22+ (a;) — ► 1 and a + {uj) — > 1 at co> — > 00, the series ( |A.ll ) is absolutely convergent in the 
domain Su; < n, and we arrive at 

[7 21+ M= ^ 22+(mK(m)e '" +Q(e- 3 "), 9*<7r. (A.12) 
u — in 
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Substituting (|A.12|) into the equation for U 22+ , putting uj = —in and shifting the integration 
contour to the lower half-plane, we obtain 



7T . 



U 2 2+{-i^) = U n+ (m) = 1 + -U 2 2+(iTr)e 
where (C « 0.57721 ... is Euler's constant) 

7 = -i d z \og a+(z)\ z= 



-2nt-2 



(t + 7 + ^)+o(e" M ), (A.13) 



-(log2-C) . 

7T 



Similarly, we analyse the equations for Un + and U12+ and obtain 



1 



C/ n+ (-i7r) = U 22 +{m) = 1 - -tW^e* 1 + o (e^) 
EW-«r) = -U 21+ (iir) = -7te-^-\t + 7) + o (e~ 2 ^) . 



Finally, using relations flO) ) with Eqs. flO^D , (|Q3|) and (|X7T5|) we find 



l+p 2 (q+l) p 
qp 1 — p 2 



+ o {p 2 ) 



where p and q are defined as 



P=- 2 e 



-vrt-1 



g = 27rt + 2(log2 - C) . 
Since the estimate ( A.12| ) is uniform in the domain < it, we have with ( |A. 1| ) 
$ 21 _(0) = U 21+ (0) = -U 2 2 + {in)e- t ^ 1 + O (e- 3rf ) . 



(A.14) 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



Comparing with ( |A.8| ) and using ( |A.16| ) we find sin<^ = 2p, and thus we arrive at (|5.9|) . 



A.2 S > 1 



The asymptotic analysis of the RHP ( }4.6| ) for S > 1 is quite similar to the case 5 = 1. Now 
the boundary values of the matrix $(u;) introduced in ( |5.13| ) on the real axis satisfy 



(A.19) 



with 



sinh 



GJu) 



c^ s \uj) sinhwS 



2_g— iut 



-a 



e 2 



sinh 



(A.20) 



sinhoj(S-l/2) 
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As for 5 = 1 this matrix can be factorized into the product of matrices M + and M_, which 
now are 

sinh — 1 sinh — 

M+{u>) = I - a? \co) — -^e-V_, M_(w) = / - , 2 V + (A.21) 

smhtUD a_ (u;) smhcui? 

The following considerations almost completely repeat the corresponding part of the sub- 
section S — 1. First, we again come to the new RHP on the contour T. The only difference is 
that now a G (0, £). Thus we obtain that $(u;) ~ 7 up to exponentially small corrections every- 
where exept w6l. In order to improve this estimate we use the corresponding singular integral 
equations. Now the asymptotic behaviour of is defined by residues of the integrands in 
the points uo = ^k. Respectively the solution of the RHP is given by the asymptotic series of 
exp{ — ^k}. Note also that just like in the case 5 = 1 the matrix $(0) can be parametrized by 
a single parameter. The explicit expression now has the form 

(cos 09 + 7? sin 09 sin 09 , 
|. (A.22) 
r] cos 09 — sin 09 COS 09 
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Figure 1: Magnetic susceptibility of the doped S — 1 (a) and 5 = 3 (b) system as a function of 
the concentration x h of holes (full line). The broken line is the contribution of the 'background' 
mode e 2 s to x- 
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Figure 2: The jump-contour for the new RHP ( A.2|) 
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Figure 3: Wilson ratio ( |6.1|) as a function of the hole concentration Xh for the S = 1 and 5* = 3 
system. 
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